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ABSTRACT 


An analytical theory of electromagnetic wave scattering from an inhomogeneous 
medium with a slightly rough boundary surface is formulated. The inhomogeneity in the 
medium is assumed to vary continuously in the vertical direction. In addition, it is 
also assumed to have a small random variation in the horizontal direction. The medium 
is assumed to consist of two layers. Maxwell's equations are solved by using the small 
perturbation method together with Fourier transform technique. The resulting differential 
equations are solved by using WKB and variation of parameter methods. Field amplitudes 
in each medium are determined by taking boundary conditions into account. The expres- 
sions for first order polarized radar backscatter cross-section ct° are obtained. 

An attempt is made to apply the developed theory to compute sea ice scatter. 

The complex permittivity of sea ice, which depends on both the temperature and 
salinity, varies with the depth of sea ice. In addition, there is certainly some variation 
in the horizontal direction. Thus, the developed model may be able to give useful 
estimates when applied to sea ice scattering. Numerical calculations are performed 
for polarized radar backscatter cross-section (oyy° and ct^ 0 ) at two frequencies, 

13.3 GHz and 400 MHz. It can be shown that WKB method is applicable at both of 
these frequencies. These theoretical results are compared with the experimental re- 
sults obtained from NASA Earth Resources Program mission 126. Theoretical results 
give the same absolute value of a° and the relative variation among the six ice types 
as is given by the experimental results. 
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1.0 INTRODUCTION 


A theoretical model has been developed to describe the backscatter from 
sea ice. Computations based on this model agree well in principle but not always 
in detail with experimental observations performed in 1970 over the Arctic Ocean, 
north of Alaska. The theoretical model presented here takes into account small- 
scale roughness of the surface, average variation in dielectric properties with 
vertical position inside the ice, and horizontal small-scale inhomogeneities within 
the ice volume. Experimental observations with which the theory is compared were 
made by NASA Earth Resources Program aircraft in April, 1970, using scatterometers 
at 0.4 GHz and 13.3 GHz as well as an imaging radar at 16.5 GHz. 

At the two higher frequencies, both theory and experiment indicate that 
the thinnest ice gives a moderate strength return. Ice of the next thickness cate- 
gory (in the neighborhood of 5 centimeters) gives the weakest return, and ice with 
increasing thickness results in increasing signal return. Both theory and experiment 
also show at 0.4 GHz that the thinnest ice gives the weakest return with the signal 
increasing as the ice thickness increases up to about 1 meter thickness and then de- 
creasing for ice thicker than 1 meter. 

This paper is a brief summary of some of the material contained in the doctoral 
dissertation at the University of Kansas by S. K. Parashar.* 


2.0 THEORY 


Ice on the polar seas is a complicated medium. The surface is exposed to 
blowing salt water and low air temperatures; the bottom is exposed to the sea. As the 
ice freezes, the salt water concentrates into brine and forms small pockets within the 
crystalline structure. Eventually, the pockets of brine also freeze if the temperature 
becomes low enough. Diffusion causes migration of the salt away from the surface to- 
ward the bottom of the ice; but some of the salt at the surface is replenished by the 
blowing sea-water spray. The amount of replenishment on the surface, the amount 
of freezing, the amount of snow cover and the rate of diffusion are all functions of the 
time of year and of the shorter-term weather variations. 


*S. K. Parashar, Investigation of Radar Discrimination of Sea Ice , The 
University of Kansas, 1974, Ph.D. Dissertation. This material was also presented at 
the URSI Specialist Meeting held 23-26 September 1974 in Bern,Switzerland . 




Because of the different” temperature regimes at the surface and bottom of the 
ice, and because of the diffusion processes that take place, the dielectric properties 
of the ice vary with depth. 

Because of the way in which the crystals are formed and of the way in which 
the brine concentrates in the ice, a vertically-oriented structure appears in cross- 
sections of the ice. Thus, one expects larger random inhomogeneities in the horizontal 
direction, but larger regular variation in the vertical direction. The surface of the 
ice is initially smooth but becomes rough because of fracturing, because of precipita- 
tion effects, because of the effects of spray hitting the surface and freezing irregularly, 
and because of the normal irregularity associated with the freezing of the brine. The 
surface effects due to the freezing precipitation and spray are relatively small scale; 
but some of the mechanical effects can be very large, as in the case of pressure ridges. 

An additional phenomenon, which is ignored in the theoretical model developed 
here, is "rafting." This results when one sheet of ice is pushed beneath another due 
to wind and current in the sea. The result is discontinuity within the ice where the 
bottom of one sheet that formerly was in contact with the sea adjoins the top of the 
other sheet that formerly was in contact with rhe air. Because of the complexities 
associated with rafting, the theoretical model ignores these effects. 

Figure 1 illustrates the fundamental theoretical model for this paper. Space is 
divided into three major regions: atmosphere, sea ice, and sea water, with the sea ice 
further subdivided into two other regions. Medium 1 is a region of decreasing permit- 
tivity, and medium 2 is a region of increasing permittivity with z = z^ the boundary 
at the level of minimum permittivity. In this model the permittivity in the ice is 
assumed to consist of the sum of two components; of these, one is a function only of 
vertical position within the ice and the other is a function only of horizontal position. 
The vertical ly- varying component has a smooth variation and the horizontal variation 
is random. The surface of the ice adjacent to z ~ Zq is considered rough, and the 
roughness at the bottom, z = , is ignored. Of course, the bottom of thick sea ice 

is known to be very rough, but the signals are not expected to penetrate with enough 
strength that this roughness needs to be considered. 

Figure 1 also illustrates that the incident angle is 0 and the scattering angle 
with the vertical 0^. The incident wave is in the XZ plane, but the scattered wave 
goes in a plane through the Z axis, making an angle 0 $ with the X axis. Similar angles 
may be defined within the sea ice. 
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Figure 1 . TheoreMcal Sea Ice Model . 
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In the analysis that follows, the profile within the ice of and is bi- 
linear; that is £q-| is a decreasing linear function and is an increasing linear 
function. The values at the top and bottom and at the central junction point, z, 
are based on temperature and salinity values for the type and thickness of ice being 
considered as reported by Hoekstra and Cappillino, 1971. Thus, the bHinear pro- 
file is an approximation to the curvilinear profile actually measured by other invest* 
igators. No such profiles for the ice were measured at the time of the backscatter 
measurements. 

Since no good information is available on the horizontal inhomogeneities 
either within the ice or at the surface, arbitrary selections of the parameters for 
these variations were made when the data were compared with the theory. 

The waves in the homogeneous media air and sea satisfy the usual equations: 

V z £ + K Z E = O 0) 


Co 1 AJLG- 


In the inhomogeneous medium, a perturbation solution results in successive 
iterations. The wave equation to be solved is 


y L e'cx,W; 


•_ E 
if/*) 


where k 1 is the wave number in the ice, E 1 is the electric field in the ice, and e‘ is 
the permittivity in the ice. The value of k 1 is given by 

K 2 Cx,y, Z) = + < 4 > 


where k f Q is the unperturbed value and ^‘(x,y) i s the randomly varying component 
associated with the first order solution. Thus the total field, after iterations, is 


E - <=, +- E^-E - 
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Considering the equations for the zero order and first order solutions, we 
observe for the homogeneous medium that both orders of solution must satisfy the 
usual homogeneous wave equations. 


v 2 e 0 + = ° ' 

? 2 E, + «?• E, = O (7 

Equation (8) shows the zero order form of the wave equation in the inhomogeneous 
ice medium, and (9) shows the equation for the first order solution in this medium. 
Clearly, (9) involves substitution of the solutions of (8) into the equation for Ej'. 


V 2 Eo -t C ?0 to — 


C £=o » Eo ) 

e'dCV 

1 1^ 2. . 


(8) 


7 z e,'+ = - E ° IS 


+ V 




-V 


Cve,c^,y)-Eo3 __ v 




e'GO 


(9) 


known t 

Ex, 0 


The solutions of the equations in the atmosphere (6 and 7) are of course well” 
For horizontal polarization, they are 



E 



-ikxSir.© f jk 2Cos© , -^KZCos© 

e Le +vc 


( 10 ) 


where V is the reflection coefficient. 


5 



In the homogeneous medium (atmosphere), we may write the solutions of 

equation (7) as a summation of plane waves in terms of their two-dimensional Fourier 

transforms A , A , and A , as shown in 
x' y' z e 

E x . C *,y, z) = ±- j[ A x e' CKxX + ^ d. K x d.k y 

71 -v (11) 


-°<3 




— oo 


(12) 


— chO 


(13) 


where 

K-, 


C K z - e-KS) 


V; 


(14) 


The zero order solution in the inhomogeneous medium can be obtained by 
solving the homogeneous wave equation 


V^E o + <,(.£) E 0 = O 
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This results in 


t 




E 2 — O 


(16) 
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K,cz) = LkoCz)-<x] , a=Es 


in Q 


(18) 


This solution was obtained by use of separation of variables, and the solution takes 

into account the vertical variation of k 1 by use of the WKB method. 

o 

The first order solution is obtained by solving an inhomogeneous wave equation 
for the Fourier transform of the field component. 


fjC d K, C 2 .) — 2Q 

dz 2 

where S is the Fourier transform of 

L^' Cx ^ )E ^] 

C>Z: £! 0 


This equation may be solved by methods similar to that used for solving the homogeneous 
equation (15) and variations of parameter methods (Boyce and Diprima, 196 7), but, of 
course, the solution is much more complicated, as indicated by its form: 


FC^X; 


JW&) 


&,Oc x ,ts y )e 
- Ca'M*' 


}L 


+ 


+ B z Ck^.,K y )e *-=> 

i f S z C^,Ky,it) 


2 3 

e. 




e 

£ 


; Zo 


i dz' 

J' J K'c^O / 

, / 

dz 


-if ( e 5>cK x ,K^z^ e. J *o 

■ ’ L- j*\&> 


7 


( 20 ) 



(21) 


Here, 


K,'ca5= [kies? 


To obtain the complete solution, the boundary conditions at the interfaces 
must be matched. The boundary conditions at the smooth interfaces are, of course, 
straight-forward. However, some discussion of the boundary conditions at the rough 
surface is in order. Since the interface between air and ice is rough, even hori- 
zontally polarized waves may have some Z component of the electric field in the re’ 
fleeted and transmitted values. Thus, the local slope of the small perturbations of 
the surface height z(x,y) must be considered. Hence, the boundary conditions that 
must be satisfied at this interface are 


V E 3 = 


"5 g O VU) rr _ F ' N 
75U fc *) 


Ev-E x = - 


Hy-VK, 


H x - H y 


J 


Here, we assume that the maximum value of z(x,y) is small enough so that 


(22) 

(23) 

(24) 

(25) 


J K 0 Cz';a,z ~K 0 (o^zc^y) 


(26) 


And, also, this means 

KoOcx.y)] ~ (27) 

The total solution is very complicated and is not presented here in its entirety. 
Rather simply, the general form is given here. Equation (28) shows the form of the 
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horizontally polarized solution and 






| Rvj CKS(n©j, 0)| A-TT K Cos 2 © 


(28) 


The vertically polarized solution is of the form 

| — A x \y) q 7 O) Cos © — A^CKSj-n©;, o) S\t70| J^> 

- lZ ^ z 'Cos z e (29) 

O . - o 

^VH - ^HV ~ ° 

Cross~poIarized values of this differential scattering coefficient are zero in this 
solution. 

The coefficients Aj are the Fourier transforms of the field solutions and be" 
cause of the assumptions in obtaining the solution, they may be expressed as the sum 
of components due to surface roughness and due to horizontal inhomogeneities for 
ths 1 permittivity in medium 1 and medium 2. That is 


ft l C<SiT 10,0) = + M 2 - e z +N vO 

Z - F.T. Of 2C*,y)y €, = F.T. Of fe,Cx,y) 

F.T. of Z 2 Cx,y) OD 


Here the constants M. include the effect of the smoothly varying permittivity change 
with depth in the ice. If we assume an exponential form for the correlation function. 
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the solutions of (28) (29) can be expressed in terms of the correlation lengths by 


<r°acsi-r,e,d)= m 0 m* r m? 

L 0 + 4 14 Siv, 2 ©)^ 

, MzM? , fiyOTv^MsM* ~j 

0 + 4^5^^ 0+4i%s^e)\|- < 32 > 

- 4-TT K 1 Cos Z e 


where L $ is a normalized surface correlation length in layer 1, and l_ 2 v is a normalized 
volume correlation length in layer 2. The standard deviations for the respective media 
are given by a z , a^, o^y' respectively. 

This outline of the theory has only attempted to give an overview. The details 
require more space than is available for this paper. Nevertheless, the general con- 
cepts of the model, based on known variations of the permittivity with depth for 
different thicknesses and types of ice with random inhomogeneities in the volume and 
on the surface, give a solution of the form shown. 


3.0 COMPARISON OF THEORY AND EXPERIMENT 


The experiment was performed in April of 1970 directly north of Point Barrow, 
Alaska, over a region of the Arctic Ocean containing a wide variety of first-year 
and multi-year ice. The P-3 aircraft contained two scatterometers, one at 0.4 GHz 
made measurements with vertical, horizontal and cross polarizations. At 13.3 GHz, 
the scatterometer measured only vertical polarization returns. At 16.5 GHz, a side- 
looking radar was used to produce images with vertical, horizontal, and cross polar- 
ization. 

The flight line was laid out as shown in Figure 2. The scatterometers were 
flown in parallel lines across an area viewed from four directions by the side-looking 
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MISSION 126 - MEASUREMENT OF SEA ICE BY RADAR 


INSTRUMENTS; 



oAOO MHz SCATTEROMETER 
013,3 GHz SCATTEROMEFER 
016,5 GHz $UR 



Figure 2 



radar. That is, the side-looking radar images were produced by pointing the radar 
toward the inside of a square flight track, and all of the scatterometer lines were 
within this square . 

The initial categories used in analysis of the scatterometer returns are listed 
in Table 1 . Categorization of the ice observed by the scatterometer was achieved 
by interpretation of stereo photographs taken at the same time as the scatterometer 
data. The techniques for analyzing the photographs were shown to us by Mr. Vern 
Anderson, of the Photographic Interpretation Corporation, Hanover, New Hampshire. 

Separation into one of the seven categories was achieved from the photographs. 
The assignment of thicknesses is based upon normal thicknesses for ice in such cate- 
gories. No direct measurement was made of the ice thickness. 

Categories 2, 3, and 4 all represent quite thin ice, since the maximum thick- 
ness of category 4 is only 30 cm. Categories 6 and 7, however, represent quite thick 
ice. In general, the boundary between ice too thin to be of significance to navigation 
and ice that will impede navigation is somewhere in the 5, 6, or 7 categories. On 
the other hand, even new ice can have a significant effect upon heat transfer between 
the water and the air. 

Figure 3 shows the 13.3 GHz experimental results. At vertical incidence 
open water gives the maximum return, although at angles beyond about 20°, it might 
be confused with some other category; and at angles of 50 or 60 degrees, it gives one 
of the weakest returns. Category 7, multi-year ice, gives a stronger return than any 
other ice type at most angles of incidence. And category 6, the thick first-year ice, 
also clearly gives a strong return. Category 5 gives a fairly weak return. Figure 4 
shows theoretical computations for the same frequency and ice conditions. The shapes 
of the curves are somewhat different from those of Figure 3. However, the order is 
almost the same at angles beyond about 25°. Absolute magnitudes predicted by the 
model are weaker at the larger angles than experiment shows. Perhaps a modification 
of the parameters for the surface roughness would result in a less steep slope and con- 
sequently larger signals at the large angles, like those observed experimentally. 
Nevertheless, the theoretical computations agree reasonably well in order with the 
experiments. 

Figure 5 shows the corresponding theoretical calculations for horizontal polar- 
ization. The order of the categories is the same as observed experimentally with 
vertical polarization, but the slope of the horizontally polarized theoretical curves 
is muchsteeper than that of either the vertically polarized experimental or theoretical 
curves. Unfortunately, no experimental data were available for horizontal polarization 
to compare with this theory. 12 



TABLE I. 


CATEGORIES USED IN SCATTEROMETER 
INITIAL ANALYSIS OF SEA ICE 


CATEGORY 

1) Open Water 

2) New Ice 

3) Thin Young Ice 

4) Thick Young Ice 

5) Thin First Year Ice 

6) Thick First Year Ice 

7) Multi-Year Ice 


THICKNESS (CM) 
0 

0-5 
5-18 
18 - 30 
30 - 90 
90- 180 
180 - 360 
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Figure 3. Experimental ct° vs . 0 for Different Categories, 
13.3 GHz, W. 


ORIGINAL PAGE IS 
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Figure 4. Theoretical a° vs. 0 for Different Categories, 
13.3 GHz, VV. 
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Figure 6 shows the experimental observations for vertical polarization at 
400 MHz. Here the water is not as distinguishable from ice as it was at 13 GHz, 
and the difference between the different ice types is not clear. The strongest 
signals at large angles are from the moderate thicknesses of ice and the weakest 
are from the thin 5 to 18 cm. region. Note also that the absolute level of the 
scattering coefficient at large angles for this frequency is of the order of 20 dB 
lower than the absolute level at 13 GHz. 

Figure 7 shows the theoretical calculations for the same frequency as Figure 
6. The absolute level of the expected return is comparable with that observed. The 
relative sizes of the returns of the different categories at large angles is also about 
as observed. On the other hand, the shape of the curves for smaller angles does 
not seem to be the same as the observation. In general, however, the agreement 
is reasonably good between theory and experiment for this frequency and polarization. 

Observations were made with other polarizations at 400 MHz, but the results 
are sufficiently similar to those for vertical polarization that they are not presented 
here. 


4.0 SUMMARY OF EXPERIMENTAL RESULTS 


A major purpose of theory and experiment on radar return from sea ice is 
determination of thickness. Thus, the best way to illustrate the usability of the 
scattering coefficient for thickness measurements is to plot scattering coefficient 
versus the best estimate we can make of the thickness. This has been done by using 
the mean thickness in each of the categories to represent that category, thus, obtain* 
ing one thickness value for each category. Figure 8 summarizes the results at 13.3 
GHz with vertical polarization. Clearly, there is an ambiguity between new ice, 
category 2, and some of the thicker ice. However, if some other means can be used 
to establish the extent of the new ice, the remaining ice categories can be dis- 
tinguished by the strength of the scattering. In fact, on an imaging radar, the 
texture of the pattern for new ice is quite different from that for the ambiguous 
thicker ice; so determination of ice thickness does appear feasible using this fre- 
quency and polarization. 
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Incidence Angle, e (Degrees) 


Figure 5. Theoretical ct° vs. 0 for Different Categories, 
13.3 GHz, HH. 



INCIDENCE ANGLE, * DEGREES 

Figure 6. Experimental cP vs. 0 for Different Categories, 
13.3 GHz, HH. 
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The vertically polarized returns at 0.4 GHz are summarized in Figure 9. 

Note that category 2, which gave a moderately strong return for 13.3 GHz, here 
gives the weakest return. The strongest signals come from ice in the neighborhood 
of 1 meter thick, and older ice gives a weaker return. This was also the observ- 
ation on the theory. Presumably, the difference in this property between 0.4 GHz 
and 13.3 GHz has to do with the relative degree of penetration of the ice in the 
two frequencies. 

Figure 10 illustrates that the behavior of the cross-polarized return is similar 
to that of the like-polarized vertical return. Behavior for horizontal polarization 
is similar and is not reported. 

Estimates were also made of the ability to distinguish ice categories on the 
16.5 GHz imagery. Because the gray scale rendition of the images was not good, 
only three categories of ice plus the open water were used in the image analysis. 

That is, categories 2 and 3 were combined, categories 4 and 5 were combined, and 
categories 6 and 7 were combined. The resulting categories are shown in Table 2. 

With this type of categorization over 80 percent "correct" identification was achieved. 
With this good agreement between the radar result and the photo- interpretative cate- 
gories, it is possible that a significant part of the disagreement occurs because of 
misclassification on the photographs rather than on the radar, although, of course, 
this cannot be quantitatively established. Figures 11 and 12 show images of the ice 
for the vertical and cross-polarized returns, with the different categories identified 
on the illustrations. 


5.0 CONCLUSIONS 


The theoretical model developed here gives results that coincide reasonably 
well in order of magnitude and in order of responses from the different types of ice. 

The theoretical responses calculated vary somewhat more rapidly with angle of in- 
cidence than do the experimental responses, with the difference being greater at 
13,3 GHz than at 0.4 GHz. 

The results of both experiment and thec-y indicate that the 13.3 and 16.5 GHz 
systems appear likely to be successful for measuring ice thickness. Presumably, the fact 
that there is moderate agreement between theory and experiment over the entire fre- 
quency range studied means that we can use the theory to interpolate for results at 
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TABLE II 


CATEGORIES USED IN IMAGE AND LATER 
SCATTEROMETER ANALYSIS OF SEA ICE 

CATEGORY 

A) Open Wafer 

B) Thin First— Year Ice 

C) Thick First— Year Ice 

D) Mulfi-Year Ice & 

Very Thick First— Year Ice 


THICKNESS (CM) 
0 

0-18 
18 - 90 
90 - 360 
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FREQUENCY • 400 MHi 
POLARIZATION • W 
SITE 93 



Figure 9. Experimental a° vs. Ice Thickness for Different 0, 
400 MHz, VV. 



Figure 10. Experimental <P vs. Ice Thickness for Different 0 
400 MHz, VH . 






other frequencies. However, the frequency spacing is so wide that interpolation 
may not be very reliable because the size of the important random variations in 
different parts of the frequency range probably also differs. 


REFERENCES 


1 . Hoekstra, P. and P. Cappillino, "Dielectric Properties of Sea Ice and Sodium 
Chloride Ice at UHF and Micro-Wave Frequencies,” Journal of Geophysical 
Research , vol . 76, no. 20, pp. 4922-4931, July, 19711 

2. Boyce, W. E. and R. C. Diprima, Elementary Differential Equations and 
Boundary Value Problems , John Wiley and Sons, !nc., New York, May^ 1967. 


22 



